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For a group of charged particles obeying quantum me- 
chanics interacting with an electromagnetic field, the 
charge and current density in a pure state of the sys- 
tem are expressed with the many-body wave function 
of the state. Using these as sources, the microscopic 
Maxwell equations can be written down for any given 
pure state of a many-body system. By employing semi- 
classical radiation theory with these sources, the micro- 
scopic Maxwell equations can be used to compute the 
strong radiation fields produced by interacting charged 
quantal particles. For a charged quantal particle, three ra- 
diation fields involve only the vector potential A. This is 
another example demonstrating the observability of vec- 
tor potential. 



Five radiation fields are perpendicular to the canonical 
momentum of a single charged particle. For a group of 
charged particles, a new type of radiation field is pre- 
dicted to be perpendicular to A (x J ; t) x [V x 
where <P is the many-body wave function. This kind of 
radiation does not exist for a single charged particle. The 
macroscopic Maxwell equations are derived from the 
corresponding microscopic equations for a pure state by 
the Russakoff-Robinson procedure which requires only 
a spatial coarse graining. Because the sources of fields 
are also the responses of a system to an external field, 
one also has to give up the temporal coarse graining of 
the current density which is often supposed to be critical 
in the kinetic approach of conductivity. 



Copyright line will be provided by the publisher 



1 Introduction This paper discusses three closely re- 
lated problems for a group of non-relativistic charged par- 
ticles obeying quantum mechanics (QM): (1) the radiation 
fields produced by the current density in a pure state of 
the many-body system; (2) derivation of the macroscopic 
Maxwell equations from the corresponding microscopic 
equations in a given pure state; and (3) the consequence 
of the averaging procedure used in (2) to computing the 
electrical conductivity. 

In condensed phases, there exist abundant radiation re- 
lated phenomena. Cherenkov radiation and transition radi- 
ation are two prominent examples! 1]. Although quantum 
electrodynamics (QED) is believed to be applicable to all 
radiation problems, it is only feasible to calculate the weak 
field (a few photons) produced by a few charged particles 
with perturbation theory lUEHUEHEl] . QED is difficult for 
the long time evolution J8[, the reaction of the atom on the 
applied field l9l [T0ll . the effect of boundary conditions iTTTIl . 



bound states and coherent radiation from many charged 
particles. To further explore radiation phenomena in con- 
densed phases, we need an alternative method for comput- 
ing field. 

The "neoclassical" or semi-classical radiation the- 
ory (SCRT) of E. T. Javnes lfl2l is concerned with the 
strong radiation of a relativistic or non-relativistic charged 
particle. The motion of the particle is treated quan- 
tum mechanically while the electromagnetic field is 
treated classicallvl8ll9lfT0lfTT1fT2l. The SCRT has suc- 
ceeded in treating spontaneous emission, absorption, 
induced emission! 13, 14], photo-dissociation, Raman 
scattering) 15), radiative level shifts lfl6l and absorption 
of radiation by diffusive electron ifTTll . Thus it is natural to 
extend SCRT to the radiation field caused by a current- 
carrying pure state of many non-relativistic charged parti- 
cles. In this paper, we merge results from the microscopic 
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response method (MRM1 H18II191 with concepts of SCRT 
to extend the latter to many-particle systems. 

Since the pioneering work of Lorentz l20l . many 
schemes have been suggested to derive the macroscopic 
Maxwell equations from the corresponding microscopic 
equations! 21 22 23 1. To obtain the macroscopic fields 
and sources, usually both temporal coarse graining and 
spatial coarse graining are taken in addition to ensemble 
average. By assuming that the motion of charged parti- 
cles obeys classical mechanics (CM), Russakoff|24| and 
Robinson 1231 (RR) proved that (i) only the spatial coarse 
graining is relevant; (ii) the spatial coarse graining is com- 
patible with the ensemble average 1 26]. Since the motion 
of conduction electrons in solids is described by QM, it 
is desirable to extend Russakoff-Robinson's procedure to 
such a situation. 

Dropping the temporal coarse graining has a pro- 
found consequence. In both microscopic and macroscopic 
Maxwell equations, current density and charge density are 
the sources of electromagnetic field. On the other hand, 
current density is also the response of a system to an ex- 
ternal field. Therefore Russakoff-Robinson's procedure 
means that the temporal coarse-grained average is irrele- 
vant to derive the irreversibility caused by a conduction 
process. 

In Sec|2] we write out the microscopic Maxwell equa- 
tion for a pure state. In Sec|3] we apply SCRT to analyze 
the radiation fields from the current density caused by one 
and many charged particles. For a quantal charged parti- 
cle, three radiation fields involve only the vector potential. 
Here we have one more example demonstrating the phys- 
ical reality of vector potential. In contrast to the classical 
radiation field, five radiation fields are perpendicular to the 
canonical momentum of particle. For a group of charged 
particles obeying quantum mechanics, one of its radiation 
fields is perpendicular to A(x^, t) x [V x (Vj<?')]> where 

is the many-body wave function, A(xj,t) is the vec- 
tor potential felt by the j th particle at time t. This kind of 
radiation does not exist for a single charged particle. In 
Sec|4] we apply the RR procedure to derive the macro- 
scopic Maxwell equations for a group of charged parti- 
cles obeying QM. We show that without temporal coarse 
grained average, the irreversibility of conduction process 
is still included in the present theory. 

2 Microscopic Maxwell equations for a pure state 

Consider a system (S) with N e electrons and AT nuclei 
in an external electromagnetic field (J 7 ) described by vec- 
tor and scalar potentials (A, (f>). Denote the coordinates of 
N e electrons as rj, ra, • • • , Tn b , the coordinates of AT nu- 
clei as Wi , W2, ■ ■ • , WjV. In the external field (A, <fi), the 
state ^'(i) of the system is determined by the many -body 
Schrodinger equation 



where H'(t) = H + Hf m (t) is the total Hamiltonian for 
S+ J-, H is the Hamiltonian of system, Hf m (i) is the field- 
matter interaction: 



N, 



Hfm{t)=Y J {^;[H^^)- Vi + V,-- A(r 3 .,t)] (2) 



2m 



e 2 A 2 (r,,t) 
2 m 



e0(r 3 -,i)} 



E 

L=l 



{ - l ^l[A(W L ,t) ■ V L + V L ■ A(W L ,t)} 



2M T 



(Z L e) 2 A 2 (W L ,t) 
2M L 



z L e<j>(yr L ,t)}. 



The time dependence of Hf m (t) comes from the ex- 
ternal field. The arguments of W are {tx,Ti,--- ,?n ; 
Wi,W 2 ,--- ) W A r;i). 

In QM, the state of the system at time t is described 
by the solution ^'(t) of Eq.(fl]), the microscopic charge 
density and current density j in state \P' are com- 
pletely determined by \P'(t). The microscopic charge den- 
sity p w> in state is /9*'(r> t) = J d,T^'*p{r)^' , where 
p(r) = £ 3 -eS(r _ r .) _ £ L Z L eS{r - W L ) is the 

charge density operator of the system, dr = dridr 1 is 
the volume element in the whole configurational space, 
dr 1 = dr2 ■ ■ ■ drjv e ciWi • • • dW^/. Carrying out the in- 
tegrals and using the antisymmetry of for exchanging 
the coordinates of two electrons, 



p W '(r,t) = N e e / dr 1 *'*'* -^Z L e / dr L W'*, 



L 



(3) 



where dr L = dr\ ■ ■ -dr Ni! dWi ■ ■ ■ dW L ^idW L+ i ■ ■ -dW^f, 
the arguments of the in the first term are (r, r 2 ■ • ■ rjv e ; 
Wi ■ • ■ Wjv), me arguments of the W in the second term 
are(n,--- ,rjv e ; Wi, • • • , W L _i, r, W L+1 ■ ■ ■ ,W A r;t). 

2.1 Microscopic Maxwell equations in a pure 
State If the system is in a state ^'(i), the microscopic 
Maxwell equations are 



V • b*' = 0, V x e*' = -db*'/dt, 



(4) 



and 



V • e*' = p*'/e , c 2 V x b*' = f'/eo + de*' /dt, (5) 

where e* and are the microscopic electric field and 
magnetic induction at (r, t) in pure state 

Applying d/dt on the first inhomogeneous equation in 
(0, and V- on the second, one has the continuity equation 
for a pure state ff'|27l: 



ihdV'/dt = H'&', 



(1) 



dp*'(r,t)/dt + V f (r,t) = 



t 



(6) 
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Eq. §6^ helps us find the microscopic current density [18. 1 9 1 
from : 

2m J 

N e e 2 



E 

L=l 



A(r,t) J dT 1 *'*^' (7) 



-£^A(r,f)/dr i S W > 

L=l L ^ 

where the arguments of the in the 1 st and 2 nd lines 
are (r, r2, • • ■ , rjv e ; Wi, • • • , Wjv - ; £), the arguments 
of the in the 3 rd and 4 th lines is (ri,-- - ,rjv e ; 
Wi, • • • , Wi_i,r, Win • • • j W^; i). The current den- 
sity ((TJ) is valid for an arbitrary gauge. In Eq.Q, the contri- 
bution to from nuclei is similar to that from electrons; 
we will not keep the last two terms. 

2.2 Current density and dipole moment density 
The a th (a = x, y, z) component of the dipole density 
operator is defined as 

^(r) = Y, er ja 6(r - r,) - £ Z L eW La 5(r - W L ). 

j=l L=l 

(8) 

We may extract the polarization (r, t) in state W from 



da(i)k' = J drP*'(r,i), 



(9) 



where d Q (i) | ^ / is the a t/l component of the induced dipole 
in state 

da(t)\ 9 > = j drj dT^fjr', a = x,y,z. (10) 

The time dependence of the induced dipole in state \P'(t) 
results from the time dependence of W (t ) . The time deriva- 
tive of the polarization can be found from dd a (t)\&i /dt: 

^d a (i)|*' = j dr^lf (r,t). (11) 
Combining Eqs. dlllOll U and integrating by parts, one has: 

»j£m->^(w?£.-*.»L) (.2, 

at 2m J or a or a 

-^-AJr,t) [ dr 1 ^' 
m J 

2 Ml J [W dr a * drj 



AT 

E 



( Z L e) 
M L 



A a (r,t) I dT L &'*&' 



where the arguments of in the I s * and 2 nd lines 
are (r, r 2 , • • • , rjv e ; Wi, • • • , "Wj\r; t), the arguments 
of the W in the 3 rd and 4 th lines is (ri, • ■ • , rjv e ; 
Wi,- ,W L -i,r,W £+ i-- - ,Wjv;t). Comparing Eq.tQJi 
with Eq.©, we find that the current density in state W is 
related to the time derivative of polarization in state W by: 

j*'(v,t) = dP *J t r > t) ,a = x,y,z. (13) 

Since both the free carriers and the bound electrons are 
included in Eq.©, the polarization defined in Eq.© 
contains the contributions from both carriers and bound 
electrons. 

2.3 Radiation field The state of S+J 7 is determined 
by the coupled equations (1114151 )1 14il. In ordinary materi- 
als, the induced motion of the charged particles by a weak 
external field is non-relativistic: the energy radiated by a 
particle is much smaller than the mechanical energy of that 
particle. Using direct product |t?') <gi e, b) to represent the 
state of S+J- is then allowed. 

The radiation fields produced by a current-carrying 
pure state &'(t) are! 



krad( x ; t) 



1 f , 3 , [^(x',t')] re t x (x-xQ 
47r e n C 3./ |x-x'| 2 



and 



(14) 
(15) 



{[gg^QU x (x-xQ} x (x-xQ 

lv_ir'|3 : 



where t' = t — |x — x'|/c. It is the time derivative 
dj* (x', t')/dt' of the current density that determines the 
radiation field. 

Denote Xq as the origin of the source distribution, x' 
an arbitrary point in a localized charge distribution, x the 
observation point. The integrands of Eqs. (1141151 ) can be ex- 
panded in a small parameter e = |x' — x |/|x — x |. Be- 
cause the radiation field represents outgoing energy, it must 
have asymptotic form proportional to |x — Xp | ~ 1 Il27ll28ll . 
The electric field and the magnetic induction in far region 
are 



h rad = -n x p w , e* ad = cn x (n x p w ), (16) 

where n = (x-x )/|x-x |, the arguments of ef ad , hf ad 
and p^ are (x, t) |29H30l . The vector field p* is given by 



9' 



p*'(x,i) 



1 



d 3 x' 



(17) 



Copyright line will be provided by the publisher 



4 



First author et al.: Short title 



E 

q=0 



v- 



)«^ +i r'(x',t- So ) 



where 9 is the angle between (x — Xq) and (x' — Xq). The 
q = term gives the dipole approximation, the q = 1 term 
gives the quadrupole and magnetic dipole approximation 
etc[|29][30). 

3 Radiation fields computed from 9j*'(x', t')/dt' 

3.1 A moving electron Let us consider the motion 
of an electron in an electromagnetic field (A, 0) and an- 
other external field U' (x) . An alkali atom in a dilute alkali 
gas is an example. Denote the single particle wave function 
of the particle as ^'(x, t), the average velocity v(i)|^/ = 
d[f d 3 xifj'*ri[>']/dt in state ip'(x,t) can be found from 
Schrodinger equation lfl4l[T8l : 



v(t) 



ih 

2m 



d 3 x(ip'Vip'* - ip'*Vip') 



The average acceleration a.(t)\^ 



— f d 3 xip'*A(x,t)iP'. (18) 
to 7 

dv(t)\^> / dt in state 



a(*)|^ = to" 1 y d 3 x^'*[eE - Vf/']^' 

+ At [ d 3 xtP'*(-ihV - eA)i' x B 
m 2 J 

J d 3 xtp'*tl>'(V x B), 



ihe 
2m^ 



(19) 



where the arguments of the fields are (x, t). The first term 
in Eq.dTPb is the acceleration caused by the electric field 
and external field U', the second term is the magnetic force. 
They are expected from CM. The third term is a quantum 
effect: an additional small component along the average 
direction of V x B. The ratio of the third term to the sec- 
ond is ~ a/L, where a is the characteristic length of wave 
function ip', L is the characteristic length scale in which 
magnetic field changes. 

For a single electron, Eq.Q reduces to: 



ihe 
2m 



- — A(x,iWV- 

TO 

Comparing Eq.(fT8l and Eq.(t20b. we have: 



(20) 



v(i)k 



- 1 J d 3 xf(x,t). (21) 



The average acceleration and the time derivative dj^ (x, t) / dt 
of current density are related by: 



a(t)| 



3„^(x,t) 



d 3 x 



dt 



(22) 



From Eqs. ( 117122] ) we can see that in the dipole approx- 
imation the quantal acceleration ai(t)\^i is the source of 
radiation. The radiation fields produced by the first two 
terms in Eq.dT9|> are already known in CM. The last term 
in Eq.(fT9l is the new feature of a quantal charged parti- 
cle. The corresponding radiation fields are perpendicular 
to (V x B), they should be detectable for a charged parti- 
cle moving in a non-uniform magnetic field. By testing the 
polarization of scattered fields, it can also be observed in 
the scattering phenomenon. 

Let us consider more general radiation fields produced 
by a quantal charged particle moving in an external field 
(0, A). Since we consider only non-relativistic motion, the 
effect of self-fields of the particle on its motion can be ne- 
glected. Taking time derivative in Eq.(l20b. one finds 



dfht,t)/dt = — EV>V* - — Vt/'(x)VV* 

TO TO 
2 3 

y (ihVip'*)+tp'*(-ihVip')] xB--^(Ax B) ip 'ip'' 



2m 



-A(V-/)- A^V'*(A-V)A 



+-^WW*) ■ V]A - V'W) • V]A} 



+ — 7 yvv 2 v'*+V' , *vvV-(v^'*)v 2 V''-(vV'')v 2 V'' ; 

4m^ 



2m 2 
i~h 
2m 



{ip'[A • V] (W>'*) ~ V>'*[A • V](W>')]} 



+ ^[(V^'*)A • w' - (VV/)A ■ w'*] 



ihe 2 , ,,_ , , 



ipr'*Vip')X7 • A, 



(23) 



where E(x, t) = - V</>(x, t) - 9A(x, t) /dt and B(x, t) = 
V x A(x, t) are the electric field and magnetic induction, 
the arguments of A, B, E and ip' are (x, t). Integrating 
Eq.(l23l over space J d 3 x, one can check sum rule (l22l us- 
ing integration by parts. However one observes that, ac- 
cording to Eqs. (114115b . the radiation field is determined by 
the time derivative of dy (x, f)/9i current density rather 
than by the average acceleration a(t). 

The radiation fields produced by the first three terms 
are perpendicular to the electric force, external force of U', 
and magnetic force respectively. They are well known in 
the classical description|27 i28ll29l . 

The remaining terms of Eq.d23b result from the quan- 
tum features of the microscopic current density in a 
pure state ip'. According to Eqs. ( 114115b . the radiation fields 
produced by the fourth term are perpendicular to (A x B), 
the radiation fields produced by the fifth term are perpen- 
dicular to A. The radiation fields produced by the sixth 
terms are perpendicular to (A • V) A. The ratio of the fourth 
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term to the second term is A/a, where A is the character- 
istic length of A. The fifth term is the same order of mag- 
nitude as the third term. The ratio of the fifth term to the 
second term is A/(h/ea), for ordinary magnetic field this 
is a small number. The sixth term is the same order of mag- 
nitude as the fourth term. 

The radiation fields produced by the seventh term are 
perpendicular to [(Vip') ■ V]A. The ratio of the seventh 
term to the second term is ~ (fr/ea 2 )/ B. The radiation 
fields produced by the eighth term are perpendicular to the 
canonical momentum Vip'. The ratio of the eighth term to 
the second term is ~1. In CM, no such a radiation field 
exist. The radiation field produced by the ninth terms are 
perpendicular to [A ■ V] (V^'). The ratio of the eighth and 
the ninth terms to the second term is A/a. 

The radiation fields produced by the tenth term of 
Eq.(l23b are also perpendicular to (Vip')A ■ Vtp *, where 
—iUV = rav + eA is the canonical momentum operator. 
In contrast to the eighth term, which does not depend on 
vector potential, the tenth term depends on A ■ Vip' . The 
ratio of the tenth term to the fifth term is ~1. Both the fifth 
and the tenth terms exist for a constant vector potential. 
The radiation fields produced by the eleventh term are per- 
pendicular to (ip Wtp *)V ■ A. The eleventh term is the 
same order of magnitude as the seventh term. 

Before the Bohm-Aharonov effect was discovered, it 
had been assumed that any (<fi , A ) is indistinguishable 
from ((f), A) if they satisfy A = A + VA and (f> = 
4> + dA/dt, where is A is an arbitrary scalar function of 
time and coordinate. The Bohm-Aharonov effect gives the 
first example that beside B = V x A, vector potential 
can be detected in another way: the phase shift # <f A ■ dx. 
The radiation fields produced by the fourth, fifth, sixth, 
ninth, tenth and eleventh terms involve either A or V • A. 
They demonstrate the observability of vector potential: one 
may detect these radiation fields outside a long solenoid for 
which B=0 but A ^0. 

Eq.d23T> is a resolution of dj* (x' ,t')/dt' . It directly 
results from the time dependent Schrodinger equation and 
Eq.d20b. One may wonder whether a regrouping the terms 
in Eq.(l23ll leads to a different explanation for the radiation 
fields. The answer is no: because Eqs. (1141151 ) are linear in 
Sj* (x' ,t')/dt', a recombination of terms only results to 
a different choice of the basis of fields. The new basis of 
field is a superposition of the old fields. 

3.2 A group of charged particles Taking the time 
derivative in Eq.(|7]i, we obtain 



df'(r,t)/dt=-V(r,t)n'(T,t) 
m 



(24) 



e Jk j dTW'*[-^V[/(r,x fe )-]Tv[/'(r,W a )] 

m •> k a 

{ J dT 1 [&'{ihV&'*)+&'*(-ihV&')]} x B(r,t) 



e 2 N t 
2m 2 



p3 : [A x B]n' + -A(r,t)(V -f) - -4[(A • V)A]n' 



2m 2 



h 2 eN e 
Am 2 



-^•[(W') • V]A(r,t)} 
J dr 1 ' VV 2 ^'* +f'*VV 2 !f' 
-(V^'^V 2 ^' - (V!?')\7 2 * r/ *] 



ihe 2 N e 
2m 2 



-^*[A( Xi ,i)-V](V^')]} 
J2 I rfr 1 [(V^)A(x J ^).V^' 

3 

-(W)A(x,,i)-v 3 >n 



ihN e e 2 
2m 2 

ihe 2 N, 



3 

■ J dT 1 ^A(x 3 -,t)x{[Vx(V/«p' 



+- 



2m 2 



[V x {Vj1r')\V'*}, 



where V r is abbreviated as V. n'(r,t) — N e J drf^'H/'* 
is the number density of electrons in state U and U' 
are the interaction between two electrons and the interac- 
tion between an electron and a nucleus. Because the con- 
tributions from nuclei are less important and do not have 
any new features for visible light, they are not included in 
Eq.(l24t. Of course the nuclear contributions are important 
for infrared radiation. 

There are 12 terms in Eq.d24"li. The first to eleventh 
terms are a many-body generalization of the corresponding 
single-particle terms in Eq.d23ll. The coherent scattering in 
CM is produced by several particles in which their posi- 
tions are correlated in the range of one wavelenth l!27ll28l . 
This feature is also inherited by the first to eleventh terms. 
We don't have sum rules for the system with more than one 
charged particles. 

The twelfth term of Eq. (l24b exits only for a system with 
more than one charged particle which obey QM, and rep- 
resents a new feature of the radiation field produced by the 
motion of many charged particles. If there is only one par- 
ticle, V x (Vj& ') = 0, Eq.d24l) is reduced to Eq.<|23}. The 
radiation fields emitted by the twelfth term of Eq.d24li are 
perpendicular to A (xj, t) x [Vx (Vj}P')], and are different 
to those emitted by the first to the eleventh terms. The ra- 
tio of the twelfth term to the magnetic force is A/a', where 
a' is the characteristic length of It should be detectable 
for the scattering of visible light by particles smaller than 
lOOnm. 
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In QED, the coupling between electromagnetic field 
and charged particles is treated as a perturbation. The ob- 
servable photons are expressed by external lines. To calcu- 
late a process with n photons, one needs a n th order tran- 
sition amplitude||2][3). This is impractical for a strong field 
(many photons) produced by a large number of interacting 
charged particles. In SCRT, Hf m is treated at zero order. 
Because no high energy photon appears in the usual state of 
a condensed phase and nanoscale system, Eqs. ( 114115124b 
are suitable for the emission and scattering of light ll 1211 131 

4 Macroscopic Maxwell equations For a nanoscale 
system or a macroscopic system, we usually cannot specify 
the initial condition of the system precisely in the sense that 
the state of the system in the future can be predicted to the 
maximum extent allowed by 0M ll36l . The results obtained 
in Sec|3] are not directly applicable. To explore radiation 
fields produced by a macroscopic or nanoscale system, we 
have to average the above results over a representative en- 
semble. In other words, we need the macroscopic Maxwell 
equations for a group of charged particles obeying QM. 

Suppose that the system is in a good thermal contact 
with a reservoir (£>) such that after a short equilibration 
time, the system reaches a steady state: the thermodynamic 
state of the system is specified by the intensive parame- 
ters of B (temperature, chemical potential etc.). The heat 
evolved is transferred to the bath 1 3 1 32 33 1: The motion 
of system has the same frequency as the external field. 

4.1 Spatial coarse grained average To describe 
the finite spatial resolution in a macroscopic measurement, 
for each microscopic quantity (r, t) defined in pure 
state ^'(t), one introducesl 24 25 1 a truncated quantity 

/oo 
d 3 r^'(R-r,t)/(r), (25) 
-oo 

where the scalar weight function /(r) satisfies two 
conditions: (i) d 3 rf(r) = 1; and (ii) F(k) = 

J^ao d 3 re~ lk r /(r) — > for k > fcrj- is solely de- 
termined by the type of problem and calculation we have 
in mind lEBI . Condition (ii) means that the spatial Fourier 
components of the field variables are irrelevant above 
a cut-off wave vector ho. is or any Cartesian 
component of e, b and j . Using the truncated quantity 

(R, t) is stricter than the simple coarse grained aver- 
age f}^ 1 j r& Q R d 3 r^ (r, t) for solids Il25l . where i?R is a 
physical infinitesimal volume around point R. 

The Maxwell equations for the truncated fields are 

V R • W = 0, V R xe* 7 = -db/dt, (26) 

and 

V R -i* 7 = ^ 7 /e , c 2 V R xW = W/to + dW/dt. 

(27) 



The arguments of e*", b*', j*' and p*" are (R, f), the spa- 
tial differential Vr is respect to the spatially coarse gained 
coordinate R. Because the spatial coarse graining d25l ) is 
a linear map, and Eqs. ( 1415b are linear, Eqs. ( 1261271 ) are also 
linear about all truncated quantities appear in them. 

If one distinguishes the free charges and bound charges 
in p*' and j*", various order electric and magnetic mo- 
ments of the bound charges also appear as the sources of 
fields [1241125 1 . Such a description appears in textbooks [27 
|2§1 . and is convenient for the crystalline metals, alloys 
(there is a clear distinction between free carriers and bound 
charges) and insulators at low frequency. In amorphous 
semiconductors, the hopping probability for localized car- 
riers is lT34l ~ 10 12 — 10 13 sec _1 . On a time scale shorter 
than 10 _12 sec, one cannot distinguish a localized carrier 
from a bound charge. Therefore for (i) metals, alloys and 
insulators at high frequency (inter-band transition); and (ii) 
plasma 1351 and amorphous semiconductors! 18 19.32 .33 1 
at any frequency, it is convenient not to distinguish the con- 
tributions from carriers and from bound electrons in p^' 
and j*"', i.e. not to make the multipole expansion. 

For the following three reasons, we do not need l24l 
I25ll27l to take an additional temporal coarse grained av- 
erage over fields and sources at any stage. First of all, 
Eqs. ( 1415b are linear in all quantities appearing in them: 
all the quantities are additive. Therefore for typical cut-off 
wave vector fco = 10 6 cm _1 and number density of parti- 
cles n = 10 22 cm~ 3 , the relative fluctuation of a truncated 
quantity is small lETl oc n -1 / 2 ^ 3 / 2 ~ 10~ 2 . The spatial av- 
erage d25l l is enough. Secondly, the time v^ 1 ^ 1 spent by 
a particle with speed v ~ 10 6 m-sec _1 traversing a distance 
Uq 1 is 10 _14 sec, is still in the range of atomic or molecular 
motions l27l . averaging over such a time period after spa- 
tial average (l25l > is pointless. Third, before the spatial av- 
erage d25l l. if one averages over a time period shorter than 
the time scale of atomic and molecular motions, one can- 
not eliminate fast fluctuation (27). On the other hand, if one 
averages over a time period longer than the time scale of 
atomic and molecular motions, one smeared or even elim- 
inates the scattering phenomenon HI . 

To obtain a macroscopic observable from the corre- 
sponding microscopic quantity in the pure state, we first 
truncate the microscopic quantity in a given pure state and 
then take an ensemble average. For a group of charged par- 
ticles interacting with an electromagnetic field, we can re- 
place the ensemble average with an average over all possi- 
ble initial pure states. Since for a given external field, the 
coupling between field and the system can be included with 
additional terms [represented by Hf m (t) in Eq.(f2]i] to the 
system Hamiltonian H, the state &'(t) of system at time t 
is completely determined by the value of \P' at a previous 
moment through Eq.(Q~|i. If we adiabatically introduce an 
external field, then the ensemble average is changed into 
an average over various initial values \P'(—oo) of state W . 
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Denote 



e(R,t) = X) w ^'(- 00 )]^'( R .*). 



(28) 



where £*" is any truncated quantity: or any Carte- 
sian component of e*", b*" and j*". W\&'( — oo)] is the 
probability that the system is initially in state >P'(—oo). 
W[&'(— oo)] depends only on the energy of oo), can 
be taken as either a canonical or a grand canonical distribu- 
tion (p678 of PD ). Then £(R, t) is the usual macroscopic 
observable (32). 

Because Eqs.( l26l2"71 > are linear, we can average the 
fields and sources over all possible initial pure states of the 
system. Then we obtain the macroscopic Maxwell equa- 
tions: 

V • B = 0, VxE = -dB/dt, (29) 



and 



V-E = p/e , c 2 V x B = j/e + d~E/dt. 



(30) 



Here we used the familiar symbols E and B to denote the 
ensemble average of e*' and b*" . The arguments of E, B, 
j and p are (R, t), V represents the operator respect to R. 
The macroscopic sources (p and j) and fields depend on the 
intensive parameters of a chosen ensemble. 

In view of facts that (i) the microscopic Eqs. d4l51 l are 
linear; and (ii) both the ensemble average and truncation 
d25l > are linear maps of the microscopic fields and sources, 
to obtain Eqs. (129l30l l. we can actually first take the ensem- 
ble average over the microscopic quantities, and then trun- 
cate the ensemble averaged quantities. 

The derivation suggested in this section differs from 
previous ones [21 ,22, 23 351 in two aspects: (i) to describe 
the electromagnetic phenomena in macroscopic media, 
temporal coarse graining is unnecessary (27). Only a spa- 
tial coarse-grained and ensemble average are needed; (ii) 
Because for a given external field, the interaction of the 
system and field can be written with additional terms 
in the Hamiltonian of svstem l37l . Then for each wave 
function which satisfies Eq.(Q]i, one may introduce the 
microscopic response for that state. Thus the ensemble 
average Eq.(f28b can be delayed to the final stage|T8 32 1 
rather than taking at the outset l37ll2l1l22ll23l . 

Since Eqs. ( 1415b have the same structure as Eqs. ( 129130b . 
to obtain the radiation fields produced by a mixed state, we 
only need to replace d'^ 1, (x', t')/dt' with <9j(x', t')/dt'. 

4.2 Polarization, current density and the source 
of radiation field For a group of charged particles obey- 
ing CM, the induced dipole moment is determined by the 
displacement of charges. The current density is propor- 
tional to the velocities of charges. The radiation field is 
caused by the accelerations of charges. We will show that 
for a group of charged particles obeying QM, there exist 
similar relations among polarization, current density and 
the source of radiation field: the current density is time 



derivative of polarization, and the source of radiation field 
is time derivative of current density. 

Eqs. dl3ll4ll5l ) have shown that the time derivative 
relation is correct for the microscopic quantities in a 
pure state. It is easy to see that the time derivative re- 
lation is also correct for the macroscopic quantities. For 
the truncated quantity P^R, t) = d?rP*' (R - 
r, i)/(r) and the macroscopic polarization P(R, t) = 
Es" W W '(-oo)]P* 7 (R, *). Eq.O implies 



dt 



(31) 

Sometimes, the second equation in (fJTJ is taken as a redef- 
inition of polarization! 35, 38 1 . Since Eqs. d29l30b have the 
same structure as Eqs. ( 1415b . the macroscopic radiation field 
is determined by dj a (R, t) /dt. 

4.3 Current density as the response of system 
to an external field In Secl4.ll we have seen that (1) if 
we average the motion of particles over the microscopic 
time scale i> -1 /sq , the averaged charge and current den- 
sity still oscillate in frequency k$v; (ii) if we average fields 
and sources over longer time scale (e.g. several periods of 
external field), the scattering phenomena are smeared{27| 
[T). Therefore we should not take temporal coarse-grained 
average even for describing the electromagnetic phenom- 
ena in a mixed state of a macroscopic system. 

On the other hand, the induced charge density p and 
the current density j are also the macroscopic responses 
of the system to an external field. Therefore both the pro- 
cedure in Sec l4.1l and the time derivative relation requires 
that temporal coarse graining should not be taken in the 
macroscopic charge density and current density either as 
the sources of fields or as the responses to external fields. 

The microscopic response method (MRM) used this 
fact as the starting point to calculate conductivity and Hall 
mobility 1 18,32, 19] . The entropy production is reflected in 
the existence of a steady state of the system, which is in 
good contact with a heat and material bath. Because the 
system is in a good thermal and chemical contact with 
a bath in an external field, after a short transient period, 
the system will eventually reach a steady state which os- 
cillates at the same frequency as the external field. In the 
steady state, the parameters characterizing the ensemble do 
not change: the evolved or absorbed energy are transferred 
to or taken from the bath. For a classical charged oscilla- 
tor undergoing forced oscillation, it is well-known that if 
the input energy is completely dissipated, the system will 
eventually reach a steady state. With perturbation theory, 
one can show that a solid in good thermal contact with a 
bath will reach a steady state in an external field ll39ll33ll . 

Temporal coarse graining is a key step in the ki- 
netic descriptions of the irreversibility [3 1,37]. For the 
processes caused by mechanical perturbations, temporal 
coarse grained average can be avoided. Faber[40| and 
Mott[41 1 conjectured an expression for the ac conductivity 
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for strong scattering based upon assumptions: (i) in the 
Kramers-Heisenberg dispersion relation for Raman scat- 
tering, put the final state and initial state the same; (ii) 
require the zero frequency limit consistent with the Green- 
wood's dc conductivity formula ll42l . this was realized by 
taking long time limit (i.e. temporal coarse graining) in the 
contracted Kramers-Heisenberg relation. If the second step 
was not used, the first step would end up with the conduc- 
tivity formula (expressed by the single electronic states) 
derived by the MRM in which no temporal coarse graining 
is taken[18, 19,32,39). For crystalline metals, such a for- 
mula was written down without derivation, cf. Eq.(13.37) 
of IS). 

5 Summary According to the semi-classical radia- 
tion theory, the radiation fields are determined by the time 
derivative of the current density. For a pure state of a many- 
body system, a strict current density expression has been 
obtained from the microscopic response method. Thus 
we have established microscopic Maxwell equations for a 
pure state. For a charged quantal particle, three radiation 
fields involve vector potential only. This is a new example 
demonstrating the observability of vector potential. Five 
radiation fields are perpendicular the canonical momentum 
of the particle. In a many-body system, one of its radiation 
fields is perpendicular to A(xj, t) x [V x (V J !f r ')], which 
does not exist for a single charged particle. We predict that 
this form of radiation can be detected. 

We have extended Russakoff-Robinson's ansatz to a 
group of charged particles obeying QM to derive macro- 
scopic Maxwell equations. Only the spatial coarse grain- 
ing is relevant. The charge density and current density 
as sources or responses is well-defined for a pure state. 
Every macroscopic quantity can be obtained from the cor- 
responding microscopic quantity by taking spatial coarse 
graining and ensemble average, cf. Eqs. ( 125128b . Temporal 
coarse graining is inadequate for the macroscopic fields, 
sources and the responses of system to field. The consis- 
tency among the induced displacement, the current and 
the acceleration requires that one should not take tem- 
poral coarse graining in the current density (transport 
coefficients). This agrees with the microscopic response 
method|18,32 ,33 19]. The entropy production is reflected 
by the existence of a steady state of the system in a good 
thermal connection with a bath. If temporal coarse grain- 
ing is taken in the transport coefficients obtained by the 
microscopic response method, one can recover the corre- 
sponding results obtained by the kinetic approaches. 
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